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The role of E+A and post-starburst galaxies 
I. Models and model results 
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ABSTRACT 

Different compositions of galaxy types in the field in comparison to galaxy clusters as de- 
scribed by the morphology-density relation in the local universe is interpreted as a result of 
transformation processes from late- to early-type galaxies. This interpretation is supported by 
the Butcher-Oemler effect. We investigate E+A galaxies as an intermediate state between late- 
type galaxies in low density environments and early-type galaxies in high density environment 
to constrain the possible transformation processes. 

For this purpose we model a grid of post-starburst galaxies by inducing a burst and/ or a 
halting of star formation on the normal evolution of spiral galaxies with our galaxy evolution 
code GALEV. From our models we find that the common E+A criteria exclude a significant 
number of post-starburst galaxies and propose that comparing their spectral energy distribu- 
tions leads to a more sufficient method to investigate post-starbust galaxies. We predict that 
a higher number of E+A galaxies in the early universe can not be ascribed solely to a higher 
number of starburst, but is a result of a lower metallicity and a higher burst strength due to 
more gas content of the galaxies in the early universe. We find that even galaxies with a normal 
evolution without a starburst have a H<5-strong phase at early galaxy ages. 
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1 INTRODUCTION 

The composition of galaxy types is different in the field (=low 
galaxy density environments) as compared to galaxy clusters 
(=high galaxy density environment). In the local universe this e.g. 
is de s cribed by the morphology-density (MD) relation JDresslerl 
1 19801 : lOemler] 1 1974T) . The higher the projected surface density 
of galaxies within an area, the higher is the fraction of early- 
type galaxies and th e lower is the fraction of late-type galaxies 
jDressleretal]|l997h . The field is rich in spirals and irregular 
galaxies, while in clusters mainly SOs, dSphs and dEs are found. 
Butcher and Oemler were the first to report a high fraction of 
blue gal axies in distant clusters as compa r ed to local ones ( BO ef - 



fect, see iButcher & Oemlen Jl984l. fl97a): lEllingson et alj diofjlh: 



iKodama & Bowerl j200lT) 7 Since galaxy clusters grow by contin 
uous accretion of the field galaxies, this difference is not caused 
by different formation histories but is due to a transformation that 
takes place during the lifetime of a galaxy l lCouch et alj 1 19941 ; 
iDressler et al.ll 19941) . This interpretation is supported by findings 
that the fraction of ellipticals remains constant in clusters of differ- 
ent redshift, while the fraction of spirals increases and the fraction 
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of S O galaxies correspo ndingly decreases with increasing redshift 
(see lFasano et alj|2000l) . When galaxies encounter a high density 
environment like a galaxy cluster the star formation (SF) is sup- 
pressed causing an evolution from late to early spectral types for 
these galaxies. This shows that environment plays a crucial role in 
the evolution of galaxies. While in earlier times it was believed that 
galaxy transformation takes only place in clusters and is related 
to the Intra Cluster Medium (ICM), today we know that galaxy 
properties like the SF activity, the Hi. content and morphologi- 
cal type depend on the local galaxy density in general and the 
ICM cannot be the only explanation to it. When spectra of the 
blue population of galaxies in clusters were investigated, most BO 
galaxy spectra showed strong emission lines typical of blue star 
forming galaxies. However some of the blue galaxies were discov- 
ered to show s pectra with strong Balme r absorption lines, but no 
emission lines l lCouch & Sharplesll 19871) . These spectra are differ- 
ent from the spectra of any normal Hubble type galaxy and were 
named E+A galaxies. In many clusters, these E+A galaxies con- 
stitute a large fraction of the blue galaxy population tearger et alj 
Il99d iPoggianti & Barbarclll996h . Some red galaxies in clusters 
were also discovered to show strong Balmer absorption lines. This 
effect is called the "spectroscopic BO-effect". X-ray observations 
have shown that E+A galaxies can be found near the edges of in- 
falling structures in clusters located in an area between early-type 
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galaxies in the central clu ster regions and late-ty pe galaxies in the 
field and cluster o utskirts jPoggianti et alj|2004) . Investigations of 
lHelmboldtl ( l2007h show that the gas content of E+ A galaxies on av- 
erage lies between gas-rich disk galaxies and typical gas poor E/SO 
galaxies. 

Several processes have been proposed as mechanisms for the 
transformation of galaxies in cluster environments. The most com- 
mon ones are mergers, harassment, gas stripping and strangula- 
tion. Harassment is an effect of tidal forces caused by a series of 
many weak galaxy-galaxy encounters with high relative velocities. 
The encountering galaxies are influenced by each other's gravita- 
tional field. They strip off stars and gas from each other's outer 
regions and weaken or even destroy their internal stabilities and 
disks. This can result in a star burst and/or the termination of SF 
on a long timescale of 1-5 Gyr. Harassment is a theoretically mo- 
tivated mec hanism, which has been observed in numerical simu- 
lations (see iMoore et al. 1 1 19981; lRichstond[l976h . It is believed to 
play an important role in rich galaxy clusters. Gas stripping de- 
scribes interactions between the gas in the disk of a galaxy and the 
hot and de nse Intra Cluster Medium (ICM) dete cted in X-ray ob- 
servations dGunn & Got j| 1 9721 : lOuilis et alfcOOOh . Some examples 
are viscous stripping, thermal evaporation, and ram pressure strip- 
ping. The effect is strongest when the relative velocities are high 
and the ICM is very dense, i.e. it is expected to take place in the cen- 
tral cluster regions. Ram pressure stripping has been observed (see 
ICavatte et al1l99d : iBravo-Alfaro et alj|200d : iGavazzi et ai]|2003h 
and has an indirect effect on the morphology of a galaxy as it leads 
to a truncation of SF on a timescale of > 0. 1 Gyr with possibly a 
small burst prior to the truncation. When SF stops in a spiral galaxy, 
its disk starts fading considerably, and the galaxy evolves towards 
higher bulge to disk luminosity fractions, hence towards an earlier 
morphological type. If instead of disk gas the low density halo gas 
gets stripped from a galaxy falling into a cluster gravitational po- 
tential, strangulation has been proposed to occur. The halo is a gas 
reservoir and feeds the SF. Its removal results in a termination of SF 
on a timescale of > 1.0 Gyr. The morphology is again influenced 
indirectly through the halting of SF. Galaxy mergers are galaxy- 
gala xy interactions which are most efficien t at low relative veloci- 
ties dToomre & Toomrelll972l ; lMihosll2004l) . Mergers are therefore 
more frequent in galaxy groups than in clusters. The relative ve- 
locities of galaxies in groups are comparable to the velocities of 
the stars within the galaxy, a condition very favorable for merging, 
while in clusters the galaxy velocity dispersion is higher by factors 
> 10. If the interacting galaxies are gas-rich, a major starburst can 
be triggered, followed by a truncation of the star formation on a 
timescale of ~ 0.1 - 0.4 Gyr. Observations about HI gas in E+A 
galaxies indicate galaxy-gal axy interactions or mergers as possibl e 
transformation mechanisms dChang eta l. 2001; Buvle et al]|2008I) . 
Enhanced merging within groups falling into clusters has been 
evolved as possible explanation of the observation that the galaxy 
population changes already at re latively large distanc es from the 
cluster center around 3 R,,„ (see IVerdugo et al1l2008t) . where the 
ICM is not yet expected to have a strong effect. While starbursts in 
the course of tidal stripping of galaxies are at most expected to be 
weak, starbursts accompanying gas-rich mergers have been found 
to be very strong in some cases. All of these transformation scenar- 
ios result in both a morphological and a spectral transformation of 
the infalling spiral-rich field galaxy population. The E+A galaxies 
seem to be a missing link in the process of spiral galaxies evolving 
into early-type galaxies and the investigation of this galaxy type 
will lead to a better understanding of galaxy evolution. With sur- 
veys like SDSS, CNOC and LCRS, E+A galaxies were also found 



in the field, which indicates that transformation of galaxies does not 
only take place in clus ters, but also in enhanced density environ - 
ments outside clusters dZabludoff etaill 19961 : iBalogh et alj|l999h . 
The respective timescales for the transformation scenarios need not 
to be the same and are still largely unexplored. The Hubble Space 
Telescope (HST) makes it possible to also study the morphological 
evolution of galaxies and compare the timescales of spectroscopic 
and morphological transformation. 

1.1 E+A Galaxies 

The spectra of E+A galaxies look like a superposition of a passive 
continuum spectrum of an elliptical galaxy (or an ensemble of K- 
type stars) without emission lines but with deep Balmer absorption 
lines typical for A stars. Therefore this class of galaxies was named 
E+A or k+a, which is only a spectral but not a morphological clas- 
sification. A lack of emission lines like in the E+A galaxies indi - 
cates a lack of actively star-forming regions dGomez et al.ll2003h . 
Usually the equivalent width of the emission line [On] at 3727 A 
is measured as an indicator of current star formation. On the other 
hand, Balmer absorption lines are an indicator for recent SF. They 
are produced by a dominance of intermediate age stars in the galaxy 
and are strongest for AO stars. A-type stars begin to dominate the 
light of a galaxy about 0.5 Gyr after SF stops and have a lifetime 
of about 3 Gyr. Their si gnature in the spectrum can best be see n 
1- 1 .5 Gyr after active SF dLeonardi & Roselll996l : |Poggiantil2004f) . 
This knowledge about the origin of the Balmer absorption lines has 
led to the conclusion that E+A galaxies are post-starforming, pas- 
sive galaxies which had their last star formati on about 1-1.5 Gyr 
ago Couch & Sharpies! fl987l : |Poggiantill2004h . A starburst, i.e. a 
short-term increase in the SF activity, preceding the halting of SF 
is required in galaxies with strong Balmer lines. The strength of the 
Balmer lines is a measure of the strength of the recent SF. 

Th e exact definitions of E+A galaxies varies among different 
authors. IVang et af]d2004) set the t hreshold for E+A galax ies at an 
EW[< Hy8, Hy, H<5 >] > 5.5 A, while lPoggianti et alj ( 120041) classify 
galaxies with EW(H5) > 3 A as E+A galaxies and propose an addi- 
tional threshold of EW(H<5) > 5 A for strong E+A galaxies. For the 
EW([0 n]), a threshold around 2-5 A is usually adopted, depending 
on the detection limit. In this work, we chose the H5-line at 4100 
A as an indicator for recent SF, because this line is less affected 
by emission filling than lower order lines. The threshold for E+A 
galaxies is set to EW(H(S) > 5 A. 

1.2 Goal of this Work 

In this work we will investigate different scenarios which are 
thought to be responsible for the transformation from one galaxy 
type to another. Our focus is on E+A galaxies, which are believed 
to be a transition stage in galaxy transformation from low density 
environment spirals to high density environment SO galaxies. For 
this purpose a grid of galaxy models is calculated with our galaxy 
evolution code GALEV, which models the photometric, spectral 
and chemical evolution of a galaxy. With starbursts and/or the halt- 
ing of the star formation we can simulate the effect of the differ- 
ent processes of galaxy transformation on the SFR. The GALEV 
code does not take into account the dynamical evolution or spatial 
resolution. We only consider here the effect of the various transfor- 
mation scenarios on the SFR and on the evolution of the integral 
spectrum and colors of a galaxy that it induces. 

Our goal is to identify possible progenitors and successors of 
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E+A galaxies and investigate their spectra, luminosities and col- 
ors, in order to understand the role of the different transformation 
processes. The paper is organized as follows. In Sect.|2]we will de- 
scribe our galaxy evolution code GALEV. Sect. [3] gives the results 
that can be drawn from our calculated grid of undisturbed and post- 
starbu rst galaxies. In Sect.[4]a conclusion is given. In a companion 
paper dFalkenberg et al.ll2009l . , hereafter Paper II) we will investi- 
gate the Spectral Energy Distributions (SEDs) of our post-starburst 
galaxy models and make a comparison with observations. 



2 MODEL DESCRIPTION 
2.1 Undisturbed Galaxies 

The GALEV code starts from a gas cloud of primordial abundance 
and a given initial mass. A given star formation rate (SFR) i/f(f) de- 
termines the total amount of stars formed in each time step. The 
mass distribution of the new stars is determined by an initial mass 
function (IMF) <p. The evolution of each star in the Hertzsprung- 
Russell diagram (HRD) is traced with a set of stellar isochrones for 
5 different metall icities. We have used the latest isochro nes from 
the Padova group ( iBertelli et al.ll994l ; lGirardi et al.l2003h that con- 
tain 26 stellar masses with 1067 evolutionary stages including the 
thermal-pulsing AGB phase. Each stage is determined by effective 
temperature r eff , luminosity Lb "' and life time. At any timestep 

bol.O 

the galaxy is described by a set of weighted isochrones. With a li- 
brary of stellar model atmosphere spectra, it is possible to calculate 
a spectrum for each isochrone. The adopted stellar library contains 
spectra from the UV to the NI R for all spectral typ es and luminos- 
ity classes for 5 metallicities l lLeieune et alJI 199a) . Summing the 
isochrone spectra for each metallicity, weighted by the SFR at the 
birth of the stars, and adding the spectra of the various metallicity 
subpopulations, we synthesize an integrated galaxy spectrum. For 
details about the GALEV code see lAnders & Fritze-v. AlvensleberJ 
J2003I) . 

For undisturbed spiral galaxy models we tie the SFR if/(t) to 
the evolving gas content G(t) with G{t) oc i//(t) and determine ef- 
ficiency parameters for spectral galaxy types Sa, Sb, Sc, Sd as to 
obtain, for a Salpeter IMF, at an age of 12 Gyr agreement with ob- 
servations in the UV, optical and NIR. The initial gas masses, and 
hence the initial value of the SFR are chosen as to obtain, after 12 
Gyr of evolution, the average observed B-band luminosities of the 
respective galaxy typ es Sa, Sb, Sc, Sd as determined for Virgo by 
ISandage et ail Jl985h . Fig. [T] shows the evolution of the SFR ij/{t) 
over time for the undisturbed Sa, Sb, Sc and Sd galaxy models. 

We stress that our Sa, Sb, Sc, Sd galaxy models are meant 
to describe spectral galaxy types. The observed one-to-one corre- 
spondence between spectral and morphological galaxy types in the 
local universe might, of course, not hold back until arbitrary high 
redshifts. 

The UV, optical and NIR colors, the spectra, the gas content 
and the c hemical abundances are in agreement with ob servations as 
shown in iBicker et alj d2004l) . However in contrast to B icker et al.1 
d2004l) . in this work the metallicity of the stellar population is held 
constant at Z = 0.008, i.e. at about half-solar metallicity. Half-solar 
metallicity is a typical value for spiral galaxies in the local and 
mode rately distant universe dTremonti et al.1 1 20041 ; IZaritskv et al.1 
1 19941) . There are several reasons why closed-box models used in the 
chemically consistent description do not appear appropriate for in- 
vestigating E+A galaxies: 1.) E+A are likely to be caused by ma jor 
mergers of two gas-rich, i.e. late-type spirals dKavirai et alj2007)) in 
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Figure 1. Time evolution of the SFR of the undisturbed Sa, Sb, Sc and Sd 
galaxy models. 

order to reproduce the large fraction of newly formed stars. Those 
mergers are by definition no closed-boxes, so that assuming closed 
box models are not justified. 2.) Other galaxy transformation pro- 
cesses in clusters are as well expected to be accompanied by mass 
loss. 3.) Bursts are likely to be associated with galactic outflows, 
also making the assumptions of closed-boxes invalid. Further sup- 
port comes from the non-d etections of many E+As in H I, implying 
low remaining gas-masses {Goto 2004; [Liu et al.|[2007h . 

The GALEV code does not take into account that galaxies 
might contain a high amount of dust, which has a reddening ef- 
fect on the integral colors. However, the majority of post-starburst 
galaxies are believed to contai n no or only little amount of dust 
dBalogh et ai]|2005l : lGotcll2004l) and therefore can be modelled ac- 
curately with GALEV. Some post-starburst galaxies, on the other 
hand seem to still contain some dust. Those galaxies, of course, 
cannot appropriately be described by our dust-free models 

Our spectra do not have enough resolution to quantitatively 
analyze any lines. While we implemented a calculation for the [O n] 
emission lines into the GALEV code, for the analysis of the Balmer 
lines we used the Lick calibration to calculate the line strength sep- 
arately (see Sect. |2.3| and |2.2| >. 

2.2 Lick Indices 

GALEV models include Lick indices, which are spectral absorption 
features in well-defined wavelength ranges. The names of the 25 
Lick indices refer to the most prominent line in the wavelength 
range of each individual index. Most indices are measured in terms 
of equivalent widths with the following definition: 

EW[A] = f 2 ™-W> dA (1) 
J a, F C (A) 

Here F\ is the flux for the index between two wavelengths Ai 
and A 2 . F c is the continuum flux, calculated from two "pseudo- 
continua" with wavelength ranges defined to the left and right of 
the central index bandpass. 

We model the Lick indices for i ndividual stars on th e ba- 
sis of polynomial fitting fu nctions of IWorthev et al.1 d 19941) and 
IWorthev & Ottavianil m9% . The index strength for the isochrones 
are derived by integrating and weighting with the IMF. With the 
isochrones from all contributions at each timestep, we get weig hted 
Lick indices for the galaxy (for detail see lLillv & Fritzel2006l) . 
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Table 1. Wavelength ranges for the Lick indices of the Balmer lines. 



Table 2. Effects on the SFR and timescales for the halting of SF for galaxy 
interactions. 



Index 


Blue continuum 


Index bandpass 


Red continuum 


W 


4827.875^1847.875 


4847.875^1876.625 


4876.625-4891.625 


Hy 


4283.500-4319.750 


4319.750^1363.500 


4367.250-4419.750 


H<5 A 


4041.600-4079.750 


4083.500^1122.250 


4128.500-4161.000 


US F 


4057.250-4088.500 


4091.000-4112.250 


4114.750-4137.250 



Interaction 


Effect on SFR 


T [Gyr] 


merger 


burst with truncation 


0.1,0.3 


harassment 


(burst and) termination 


1.0 


gas stripping 


(burst and) truncation 


0.1,0.3 


strangulation 


termination 


1.0 



In this work, we only investigate the Balmer line indices 
Hy, H(5p and in particular US a- In Tab. Q] the bandpass definitions 
of the index band passes and the pseudo-continua are listed (see 
iTrageretal .11 19981) . 

2.3 [O n] -Emission Lines at 3727 A 

We implemented into the code a calculation of the EW([0 n]) for 
the [O n]-line at 3727 A, which is an indicator of current SFR and 
is used for the definition of E+A galaxies. The EW([On]) can be 
modelled from th e SFR and B-band luminosity with the following 
relations given by |KennicuttUl992bh : 

SFR(M yr')^5-l(r 41 L([C>ii]) (2) 
and 

L([On]) ~ 1.4 ±0.3- 10 M -^-EW([On])ergs- 1 (3) 

-L B (G) 



Table 3. Comparison of the burst strengths obtained for three different def- 
initions in the literature. 



Type 


'burst [Gyr] 


r[Gyr] 


b [%] 


b s 


kt, 


Sa 


11 


0.1 


30 


0.01 


7.13 






1.0 


30 


0.01 


0.71 






1.0 


70 


0.02 


1.66 




9 


1.0 


30 


0.01 


0.60 






1.0 


70 


0.03 


1.40 




6 


1.0 


30 


0.04 


0.71 






1.0 


70 


0.10 


1.66 


Sd 


11 


1.0 


30 


0.24 


2.22 






1.0 


70 


0.56 


5.17 




9 


1.0 


30 


0.36 


2.71 






1.0 


70 


0.83 


6.33 




6 


1.0 


30 


0.68 


3.46 






1.0 


70 


1.58 


8.08 



2.4 Galaxy Transformation Scenarios 

The interactions of galaxies, as described in Sect.Q]have different 
characteristic influences on the SFR of a galaxy. Only the effects 
on the SFR will be the subject of this investigation. The dynam- 
ical and morphological transformation, which accompany galaxy 
transformation processes are beyond the scope of this work. 

The various mechanisms are realized with GALEV by induc- 
ing a sudden change in the SFR of a normal undisturbed galaxy. 
Three different ways of changing the SFR can describe the four 
important mechanisms of galaxy interaction. 1.) A truncation of 
the SFR on a short timescale, 2.) a termination of SF on a long 
timescale, and 3.) a burst, which is a rapid short-term increase in 
SFR before a truncation or termination of the SFR. Our assump- 
tion of a sudden increase is compatible with the findings from 
iBekki et"aT] d2005l) . who used a 2-D dynamical simulation of spiral 
mergers and also found a very sharp rise in SFR at the onset of the 
burst followed by an exponentially declining SFR. Possible scenar- 
ios for the decline in SFR after a burst are exhaustion of the fuel 
for SF, the heating and stirring of the ISM during the burst or the 
expansion of the remaining gas by a large scale galactic wind. We 
chose spiral type galaxies as most likely progenitors for E+A galax- 
ies. Bursts of SF at various times t bmsl , with various burst strengths 
b and decline timescales t are induced in the different spiral type 
models to simulate the different burst and truncation scenarios. Tab. 
[2] shows how the effects of the different transformation processes 
are described by our models in their impact on the SFR. 

To model a burst, the SFR of an undisturbed spiral is set to a 
value i/fburst at time f burst . The decline after the burst is described by 
an exponential law with a decline timescale t. 

lAburst 



W > fburst) 



(4) 



g'-'burst/r 

The maximum burst SFR lAbm-st together with t determines the 



burst strength b. In our models b is primarily defined as the fraction 
of the remaining gas at fburst that is consumed in the burst 



AG 

G(?burst) 



(5) 



G(fburst) is the g as mass at the beginning of the burst, AG the 
fraction of G(fburst) that is consumed for SF during the burst. Other 
definitions of burst strength also found in the literature are based 
on the increase of stellar mass S {6]l or the SFR enhancement Q. 



S 



(6) 



(7) 



with if/j as the SFR just before the burst. In Tab. [3] the burst 
strength values from these three definitions are compared for a few 
examples. 

A truncation or a termination of SFR without a preceding star- 
burst is implemented with a similar exponential law starting from 
the SF at the beginning of the truncation/termination t tmnc . 

We have calculated a grid for the spiral galaxy types Sa, Sb, 
Sc and Sd each with burst strengths of b=0%, 30%, 50% and 70%, 
onsets of the bursts at galaxy ages of 3 Gyr, 6 Gyr, 9 Gyr and 1 1 
Gyr and decline times of 0.1 Gyr, 0.3 Gyr and 1.0 Gyr. A burst 
strength b=0% stands for pure truncation or termination of SF de- 
pending on the decline time. In the following we will refer to galax- 
ies with starbursts, SF truncation and SF termination together as 
(post-)starburst galaxies, since one can describe a truncation or ter- 
mination as starburst with a burst strength of 0%. 
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Figure 2. Time evolution of the SFR of Sc galaxy models with bursts. The 
bursts beginning at 6 Gyr have different burst strengths of 30%, 50% and 
70%, but the same decline times of 1.0 Gyr. The bursts beginning at 9 Gyr 
have different decline times of 0.1, 0.3 and 1.0 Gyr, but the same burst 
strength of 30%. 
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Figure 3. Time evolution of the SFR of Sc type spiral galaxies with a halting 
of SFR on different timescales of 0.1, 0.3 and 1.0 Gyr, but with the same 
onset at 6 Gyr. 



2.5 The SFRs of Starburst Galaxy Models 

In Fig.[2]the evolution of the SFR tf/(t) for an Sc galaxy model with 
various burst scenarios is shown in contrast to an undisturbed Sc 
galaxy. 

Three bursts beginning at 6 Gyr with the same decline times 
of r=1.0 Gyr but different burst strengths of b=30%, 50% and 70% 
are compared. From the three burst scenarios beginning at 9 Gyr 
on the right, it can be seen that at a given gas consumption rate the 
decline time, i.e. the transformation scenario, influences the maxi- 
mum of the SFR in the burst significantly. At fixed r a burst with 
higher gas consumption has a larger t^bm-st- While a long decline 
time typical for mild harassment or strangulation results only in a 
small maximum burst SFR if/i, UISt , a medium or short decline time 
caused by gas stripping or by a merger produces much higher burst 
SFRs when consuming the same amount of gas. It can also be seen 
that the burst SFR significantly depends on the onset time of the 
burst. A later beginning implies that the galaxy has already con- 
sumed more of its gas content. Therefore a burst at a galaxy age of 
9 or 1 1 Gyr will not have as high a SFR as a burst at 6 Gyr. 

The scenarios in Fig. [3] again compared to an undisturbed Sc 
model, show the truncation/termination of the SF of an Sc model at 
an age of 6 Gyr on different timescales t. The models with r=0. 1 
Gyr and r=0.3 Gyr simulate a truncation after a merger or the strip- 
ping of gas from the disk. The model with r=1.0 Gyr is a termina- 
tion of SF as caused by galaxy harassment or strangulation. 



3 RESULTS 

3.1 Lick Index H<5 of Undisturbed Galaxies 

In this Sect, we investigate the evolution of the EW(H<5) for our 
undisturbed galaxy models for spectral types Sa, Sb, Sc to later 
compare them with the transformation models. In Fig. [4] we plot 
the EW(H<5) as a function of galaxy age. The horizontal line in Fig. 
[4] indicates the 5 A threshold, separating H<5-strong from non-H<5- 
strong galaxy models. Galaxies lying above this line during any 
time of their evolution would be observationally classified as HS- 
strong galaxies and therefore fulfill one of the two criteria for the 
definition of E+A galaxies. 
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Figure 4. Time evolution of the EW(H(5) of undisturbed spiral galaxies. The 
horizontal line marks the 5 A threshold. 

Fig.|4]shows that the EW(H<5) grows fast after the onset of SF 
and reaches a maximum value around 1 Gyr. The EW(H(5) of 5.5-6 
A of the different galaxy types declines on different timescales. 

It is to note that all undisturbed normal spiral galaxy models 
lie above the threshold and therefore have an H(5-strong phase early 
in their evolution lasting until ages of 3 to almost 6 Gyr, i.e. to 
redshifts z ~2-l for Sa and Sd types, respectively. The values for 
EW([On]) and for EW(H<5) for an Sa and an Sd model at 1.5 Gyr 
are listed in Tab. [4] The values for the EW([On]) are too high for 
the undisturbed galaxies in the H<5-strong phase to classify them as 
E+A galaxies. 

The values for EW(H<5) and EW([On]) between ages of 1.5 

Table 4. EW([On]) and EW(Hc5) of Sa and Sd model at ages of 1.5 Gyr, 4 
Gyr (Sa) and 6 Gyr (Sd), respectively. 



Type 


Age [Gyr] 


EW([On])[A] 


EW(H<5)[A] 


Sa 


1.5 


37.25 


5.83 




4 


26.98 


4.88 


Sd 


1.5 


42.06 


5.64 




6 


33.63 


4.97 
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and 6 Gyr are a striking result, because it indicates that galaxies 
which only follow the normal evolution, all go through a phase in 
which they have strong Hi5-lines without having a starburst. Not 
taking this result into account in analyses of high redshift galaxies 
leads to an underestimate of their emission lines, which are influ- 
enced by the Balmer absorption lines, for example the Ha and Hf}, 
and of the SFR derived from those. 

From Fig.|4]it is also obvious that a lower threshold of 3 A for 
EW(H(5) as sometimes used in literature is not appropriate because 
the undisturbed galaxy models have values of EW(H(5) > 3 A for at 
least 9 Gyr, i.e. until z ~0.4. 

3.2 Lick Index US for Post-Starburst Galaxies 

For the discussion of the E+A galaxies, we start with the investiga- 
tion of EW(H(5) and EW([On]) of our starburst models, because 
the presence of strong Balmer absorption lines and the absence 
of emission lines are the two criteria to distinguish between E+A 
and non-E+A galaxy models. In Sect. 13. 21 we investigate which of 
our models are Hc5-strong by looking at the EW(H<5) as a func- 
tion of time and in Sect. 13.41 we investigate the time evolution of 
EW([On]). 

In Figs.[5]and[6]fhe evolution of the EW(H(5) of different burst 
and truncation models are represented. Fig. [5] shows an Sd and an 
Sa model with a burst at a galaxy age of 6 Gyr as well as an Sa 
model with a burst at a galaxy age of 9 Gyr and 1 1 Gyr. For all 
four models the burst strength b amounts to 70% and the decline 
time is t=0.1 Gyr. For the two models with a burst at 6 Gyr, the 
EW(H(5) drops strongly at the onset of the burst before it grows and 
crosses the threshold of 5 A. This drop is caused by the dominance 
of young burst stars, which do not show strong Balmer lines. The 
Hci-strong phase starts 0.2 Gyr after the onset of the burst and lasts 
for 0.5 Gyr (Sa) to 1.1 Gyr (Sd). After the H(5-strong phase the 
EW(H(5) declines very fast and becomes smaller than in the undis- 
turbed galaxies. In contrast to the undisturbed galaxies, the decline 
of the EW(Hi5) is similar for the different galaxy models. The Sd 
galaxy shows the strongest H(5-line, because an Sd galaxy at an age 
of 6 Gyr has a higher gas content than the other spiral models which 
consume their gas faster and hence allows for the highest burst SFR. 
Also at an age of 6 Gyr the underlying stellar population before the 
burst has already a higher EW(H<5) than earlier spiral models. 

Comparing the three Sa models with the different onset times 
of 6 Gyr, 9 Gyr and 11 Gyr, Fig. [5] shows that the burst strength 
reflected in the strength of the H<S-lines depends strongly on the 
gas content and its evolution over time. For early bursts, the gas 
reservoir in the Sa model is still high enough to allow for high burst 
SFRs and to drive the galaxy to high values of EW(H<5) after the 
burst. At a galaxy age of <9 Gyr, the gas content of the Sa galaxies 
is very small because most of the gas has already been converted 
into stars. Due to the small amount of gas left in the galaxies, only 
a small starburst can occur at these redshifts and the models do not 
reach the HiS-strong regime any more at low redshifts. This helps 
us to understand why relatively few E+A galaxies are observed at 
low redshifts. 

Fig.[6]compares burst and truncation scenarios and bursts with 
different decline timescales. In comparison to the burst scenar- 
ios, the EW(H<5) of the truncation and termination scenarios in- 
creases immediately, without dropping first. However, the strength 
of EW(H<5) is much weaker than for galaxies with a preceding burst. 
In fact only the very strongest cases of the truncation scenario, i.e. 
the sudden halting of the relatively high SFR in Sc and Sd models, 
result in EW(H(5) values slightly above 5 A for 0.3 Gyr (Sc) and 0.5 



CO 




A Sd, 6 Gyr — 

\ Sa, 6 Gyr 

\ Sa, 9 Gyr 

\ Sa, 11 Gyr - 


6 






$4 






to 

I 

s 

m 2 











-9 





2 2 4 6 8 10 12 14 



Age [Gyr] 

Figure 5. Time evolution of EW(H(5) in galaxies with different burst sce- 
narios (Sa, burst of 70% at 6 Gyr, decline time 0. 1 Gyr; Sd, burst of 70% 
at 6 Gyr, decline time 0. 1 Gyr; Sa, burst of 70% at 1 1 Gyr, decline time 0. 1 
Gyr). 
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Figure 6. Time evolution of EW(H(5) in galaxies with different burst sce- 
narios (Sd, burst of 70% at 6 Gyr, decline time 0. 1 Gyr; Sd, burst of 70% at 
6 Gyr, decline time 1.0 Gyr; Sd, truncation at 6 Gyr, decline time 0.1 Gyr). 



Gyr (Sd), while no termination scenario drives the EW(H(5) values 
above the threshold. Compared to the burst scenarios, the Hc5-strong 
phase of the truncation scenarios starts earlier, e.g. 20 Myr after the 
onset of SF truncation in an Sd model. Fig. [6] also illustrates on the 
example of an Sd model that a longer decline time results 1.) in 
weaker Balmer lines and 2.) possibly in a longer H(5-strong phase 
if enough gas is available and therefore the burst SFR can be strong 
enough for the model galaxy to reach the 5 A threshold in the post- 
starburst phase. 

The H<5-strong phase can last up to 2 Gyr in cases with long 
burst timescales. 

The EW(H(5) strength does not simply depend on the burst 
strength, but on a complex interplay of burst strength, decline 
timescale, onset of the burst and progenitor galaxy type. Depending 
on the individual model either one of these factors can dominate the 
influence on the EW(H<5). 

To study the possible impact of metallicity effects on the 116- 
strength, we run Sd type constant SFR models with different metal- 
licities 0.0004 < Z < 0.05. The scenarios in Fig. [7] show three 
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Figure 7. Evolution of H<5 with time for Sd models with three different 
metallicities (Sd, 70% bursts at 6 Gyr, decline times 0.1 Gyr). 



bursts at 6 Gyr, a burst strength of b=70%, a decline time t=0.1 
Gyr, but with different metallicities. 

It is obvious from Fig. UJ that the EW(H6) depends on 
the metallicity of the galaxy. While the peak values achieved 
for EW(H<5) are somewhat higher for high metallicity galaxies 
(EW(FM) ~ 9 A for Z=0.05 versus EW(Hd) ~ 8 A for Z=0.0004), a 
low metallicity galaxy remains much longer in the HS-strong phase. 
At low metallicities the stars are generally hotter and brighter. 
Therefore stars of lower mass populate regions in the HRD that are 
described by A and early F spectra: For example at solar metallic- 
ity the 1.4 Mq star is the lowest mass reaching into the early F-star 
range. At low metallicity, longer-living stars down to 0.9 Mq reach 
this range and live there for about 4 Gyr. Stars with 0.9 Mq are 
much more numerous than 1.4 Mq stars as a consequence of the 
stellar IMF. The A-star temperature range in the HRD is reached 
at solar metallicities by stars > 1.6 Mq which live there for < 1.5 
Gyr. At low metallicities (Z<0.02) these temperatures are already 
reached by stars ^ 1 Mq which live there for > 3.8 Gyr. The life 
time of an A star increases therefore by a factor > 2 at lower metal- 
licities, while additionally there is a higher number of A stars, both 
acting together to cause a longer H<5-strong phase. These two ef- 
fects apparently strongly overcompensate the fact that the lifetime 
of 1.5-3 Mq stars is longer at higher metallicities. The lifetime of 
stars in this mass range at a metallicity of Z=0.05 is only a factor ~ 
2 higher than their lifetime at Z=0.0004 and these high mass stars 
are much less numerous. 

Because galaxies in the early Universe tend to have lower 
metallicities than galaxies in the local Universe and since a lower 
metallicity implies a longer H<5-strong phase, the probability to ob- 
serve a galaxy in its E+A phase increases towards higher redshifts. 
Therefore a higher number of E+A galaxies cannot be ascribed 
solely to a higher number of starburst galaxies but also to the lower 
metallicities and therefore a longer H<5-strong phase of galaxies in 
the early Universe. Higher maximal burst strengths enabled by the 
higher gas content add to this effect. 

In the GALEV code, two different Lick indices for H<S are im- 
plemented, H6 A and HS F . The index H<5 A , which is mainly used in 
this wo rk has a wider index b andpass and wider pseudo-continuum 
ranges dLillv & Fritz3l2006h . The values of EW(H(5 f ) are signifi- 
cantly lower at all times than the values of EW(H<5a). Therefore it 
is very important to check for the consistency of the EW(H<5) defi- 
nitions before comparing with observations of all the Balmer lines 



EW(H/3), EW(H(5) and EW(Hy). The EW(HcS) shows the highest 
values and is therefore easier to measure than the other Balmer 
lines. Moreover, in contrast to lower order Balmer lines the H<5- 
line is less affected by any potential emission line contributions in 
case of low level residual SF in the galaxy. 



3.3 H<5-Strong Galaxies 

In this Sect, we will explore those scenarios from our grid which 
become H<5-strong and identify the progenitors of the Hci-strong 
galaxies. 

As shown in Sect. 13.11 all burst scenarios with a burst at 3 
Gyr are H<5-strong since already the undisturbed galaxy models lie 
above the 5 A threshold at this time of their evolution (see Fig. 
|4](. For a certain type of progenitor spiral galaxy, the essential 
parameter to cause an EW(HS) above 5 A is the peak burst SFR 
'/'burst and not the burst strength. For different spiral progenitors, 
different peak burst SFRs are required in order to make them HS- 
strong after the burst. For example, all Sa galaxies with a burst SFR 
^burst <i 15—^ become H<5-strong, while for Sa models with a burst 
SFR below this value, the EW(Hr5) does not reach the 5 A threshold. 
There are a few exceptions of galaxies which have a value for (iwt 
above the threshold value for the corresponding galaxy type but do 
not become Hci-strong. These exceptions are due to the EW(H<5) of 
the underlying undisturbed galaxy models. Bursts at a later stage in 
the life of the galaxy cannot cause as strong an H5-line, not even for 
higher SFRs than bursts at earlier times, because the EW(H<5) of the 
underlying undisturbed model is already low. The later the bursts 
start the lower are the peak values for the EW(H<5). The EW(Hcf) of 
the underlying galaxy is also responsible for the fact that the dif- 
ferent galaxy types have to be examined separately and no global 
minimal value can be given for a burst SFR to cause an H<5-strong 
phase. 

Without exception all Sd models with a burst become H<5- 
strong. The maximum value for EW(H<5) of 8.5 A in our sample 
is reached by the Sd model with a burst of strength b=70% at 6 Gyr 
and a decline time r=0.1 Gyr. Almost all Sc burst models become 
H(5-strong. Only Sc burst models with long decline time which have 
moderate to low burst strength and the burst at late galaxy ages do 
not reach the threshold. Sb burst models with a long decline time 
do not become H(5-strong, while all Sb models with a short decline 
time have a H<5-strong phase. For Sb burst models with a moderate 
decline time of 0.3 Gyr it delicately depends on the beginning and 
the strength of the burst. Most Sa burst models do not become 116- 
strong. Only very few models with early bursts which have short 
decline times have a H(5-strong phase. 

For the models with SF truncation and termination without a 
preceding burst only Sc and Sd galaxies can possibly become 116- 
strong. No Sa or Sb truncation or termination model at any galaxy 
age or with any decline time reaches EW(H<5)> 5 A. Sc trunca- 
tion and termination models without a preceding burst only be- 
come Hc5-strong, when SF is truncated at early ages and on short 
timescales, i.e. harassment and strangulation scenarios without a 
burst can be excluded as possible progenitors for E+A galaxies. 
In contrast, the Sd models with short and moderate decline times 
become H(5-strong. Even an Sd model with a SF termination on a 
long timescale at early ages has an Hc5-strong phase. Tab. [5] shows 
a compilation of the H5-strong models with starbursts at 9 Gyr. 

It can be seen that the later the progenitor galaxy type, the 
more models become Hc5-strong. 
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Table 5. Compilation of H(5-strong galaxy models with a burst beginning at 
9Gyr. 



Type 


b[%] 


t [Gyr] 


Sa 


Sb 


30 


0.1 




50 


0.1 




70 


0.1,0.3 


Sc 


30 


0.1,0.3 




50 


0.1,0.3 




70 


0.1,0.3, 1.0 


Sd 





0.1,0.3 




30 


0.1,0.3, 1.0 




50 


0.1,0.3, 1.0 




70 


0.1,0.3, 1.0 



150 



50 







I 


Sd, b=70%, 6 Gyr, x=0 


Gyr 




Sc, b=70%, 1 1 Gyr, t='< 


.3 Gyr 
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Figure 8. Time evolution of EW([0 n]) of an Sd galaxy model (Sd, burst 
70% at 6 Gyr, decline time 0. 1 Gyr) and an Sc galaxy model (Sc, burst 70% 
at 1 1 Gyr, decline time 0.3 Gyr). 



3.4 [O ii]-Lines of Hr5-Strong Galaxies 

The second important criterion for a galaxy to be classified as E+A 
galaxy is to ensure that the galaxy has no current SF by measur- 
ing the EW([On]) at 3727 A. For this purpose we modelled the 
EW([On]) on the basis of the SFR and B-band luminosity using a 
relation given by iKennicut t (1992b) as introduced in Sect. 12.31 In 
Fig. [8] the time evolution of the EW([On]) of an Sd galaxy with a 
burst of 70%, an onset at 6 Gyr and a decline time of 0. 1 Gyr, and 
of an Sc galaxy with a burst of 70%, an onset at 1 1 Gyr and a de- 
cline time of 0.3 Gyr is shown. The EW([On]) increases instantly 
at the onset of the burst and afterwards decreases rapidly. 

It is obvious that the EW([On]) reflects the time evolution 
of the SFR. It is to note, however, that the EW([On]) remains at 
high values for some time after the burst and therefore extends into 
parts of the H<5-strong phase, as can be seen in Fig. [9] for the Sd 
and Sc galaxy models described above. In this Fig., the time evo- 
lution of EW([On]) is plotted together with the time evolution of 
EW(H(5). The peaks of the EW([On]) are not shown in this figure, 
because the values are much higher than EW(Hr5) during the burst 
and would make it impossible to see the EW(H(5) in the same plot. 
The horizontal line marks the thresholds for EW(Hr5) as well as for 
EW([On]). The threshold for EW([On]) is chosen to be 5 A, i.e 
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Figure 9. Time evolution of EW([0 11]) and EW(H<S) of an Sd galaxy model 
(Sd, burst 70% at 6 Gyr, decline time 0.1 Gyr) and an Sc galaxy model (Sc, 
burst 70% at 6 Gyr with decline time 0.3 Gyr). 

at the upper limit for this threshold from lDressler et al. Nl999h .We 
chose this upper limit taking into account that observed spectra al- 
ways have noise, which tends to bias observed EW([On]) towards 
lower values, while we calculated the EW([On]) from the SFR. 

The Sd galaxy in Fig. [9] is the model with the longest E+A 
phase of 0.71 Gyr in our grid. It is the same galaxy that has the 
strongest peak value for EW(H(5) in our sample. The beginning of 
the E+A phase is given by the drop of EW([0 11]) below 5 A, while 
the end of the E+A phase is defined by the drop of EW(H5) be- 
low 5 A. The duration of the E+A phase thus amounts only 64% 
of the H<5-strong phase with a duration of 1.11 Gyr. This example 
shows that the E+A phase only accounts a dramatically short frac- 
tion of the post-starburst phase. For many galaxy models, the phase 
of strong [O n] emission is even longer than the Hc5-strong phase as 
demonstrated in the left of Fig. [9] The Hr5-strong phase of the Sc 
galaxy model is already over when the EW([0 11]) drops below 5 A 
and this model therefore would not be classified as an E+A galaxy. 

It is striking that the models can clearly be separated accord- 
ing to the different decline times from looking at the strength of 
EW([0 11]). All HS-strong galaxy models with SFR decline times of 
0. 1 Gyr have an E+A phase, i. e. sufficient low [On ] -lines during 
their HS-strong phase, while all models with SFR decline times of 
1.0 Gyr and 0.3 Gyr never go through an E+A phase since they still 
have [On] emission during their HS-strong phase. 

Therefore we state that the conventional definition of E+A 
galaxies based on the absence of [0 11] emission above some thresh- 
old is too narrow to encompass the full range of post-starburst 
galaxies. In particular, galaxies with a long decline time are ex- 
cluded since they exhibit significant [On] emission in addition to 
their strong Balmer absorption lines. Possible processes leading to 
E+A galaxies can only be galaxy mergers and gas stripping (with a 
starburst), while harassment and strangulation can be excluded due 
to their long decline timescales. 

3.5 Spectral evolution of E+A Galaxies 

Figs. [TU] to [J_3] show spectra of the galaxy models with the high- 
est and the lowest peak value for EW(H(5). The highest value is 
reached by the Sd model with a burst strength b=70%, an onset at 
6 Gyr and a decline time of t=0.1 Gyr (see Sect. 13.2b , shown on 
the top of Figs.[T0]to[T3] The bottom of Figs.Q7Jto[T3]shows the Sa 



The role ofE+A and post-starburst galaxies - 1. Models and model results 9 




0.35 p 

0.3 - 
0.25 - 

S> 0.2 - 

E 

o 

CD 

j5 0.15 - 

X 

^ 0.1 - 
0.05 - , 

Hr 



2000 



4000 6000 
MA] 



8000 



10000 



Sd 1 Gyr before burst 



2000 



4000 6000 
MA] 



8000 



10000 



Figure 10. Spectra of an Sa model (burst 50%) and an Sd model (burst 
70%) 1 Gyr before the bursts (at 6 Gyr, decline times 0.1 Gyr). 
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Figure 11. Spectra of an Sa model (burst 50%) and an Sd model (burst 
70%) during the burst (at 6 Gyr, decline times 0.1 Gyr). 



model with a burst strength b=50%, an onset at 6 Gyr and a decline 
time t=0.1 Gyr, which has the lowest value for EW(H<S) of all our 
galaxy models classified as E+A. Note the different flux scales in 
all graphs of Figs.[TOJto[T3] The different stages of the galaxies in 
these figures are given in relation to the onset of the burst. The two 
graphs in Fig. [TO] show the spectra of the Sd and Sa burst models 
at an age of 5 Gyr, i.e. 1 Gyr before the burst starts. Therefore the 
galaxies have spectra of normal undisturbed Hubble type Sd and 
Sa galaxies at that age. The graphs in Fig.[TT]show the same galax- 
ies during their bursts at an age of 6 Gyr. The spectra show much 
stronger emission lines and a higher UV flux than the undisturbed 
galaxies. 

In Fig. [12] the H<5-strong phase of the two models is shown. 
The spectra show the galaxies at the age, when the EW(H<5) has 
reached its maximum value (i.e. at 6.6 Gyr for the Sd, and at 6.4 
Gyr for the Sa model) and at an age of 7 Gyr, i.e. 1 Gyr after the 
burst. In Tab. [6] the corresponding line strengths for EW(H<5) and 
EW([On]) are listed for these spectra. We recall that our spectra 
do not have enough resolution to measure the line strength. The 
values for the line strengths of EW(H<5) and EW([On]) listed here 
are calculated as described in Sect. |2.3| and |2.2| 

The spectra at 6.6 Gyr and 6.4 Gyr, respectively, still have 
high UV fluxes and the overall continuum does not yet look like 
a typical passive galaxy spectrum but instead resembles an A-star 
spectrum. In comparison, the two spectra look rather different, es- 
pecially when looking at the regions around the Lyman-a-line at 
1216 A and the 4000 A break. The Balmer lines are stronger in the 



Table 6. Values of EW(H<S) and EW([On]) for an Sd model (top) with a 
burst strength of 70% and an Sa model (bottom) with burst strength of 50%, 
both with an onset at 6 Gyr and a decline timescale of 0. 1 Gyr. 



Sd Age 


EW(Hr5)[A] 


EW([On])[A] 


6.6 


8.5 


4.9 


7.0 


7.1 


0.2 


9.0 


1.3 






Sa Age 


EW(Hr5)[A] 


EW([On])[A] 


6.4 


5.6 


11.7 


7.0 


3.4 


0.1 


9.0 


0.3 





Sd galaxy and no significant emission lines can be seen. In fact, the 
Sd galaxy at 6.6 Gyr has an EW([On]) lower than 5 A and would 
be classified as an E+A galaxy. However, its spectrum at the bot- 
tom of Fig.[l2]does not look like a typical E+A spectrum, i.e. like 
a spectrum of a passive galaxy with significant A-star features, par- 
ticularly when looking at the UV flux. This confusion is due to the 
fact that in the literature, spectra of E+A galaxies do not reach into 
the UV region, which, however, is the wavelength range with the 
strongest changes during the post-starburst phase. 

The Sa model has strong [O n] emission lines for a longer time 
after the burst, while its values for EW(H($) are weaker and drop 
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Figure 12. Spectra of an Sa model (burst 50%) and an Sd model (burst 
70%) in the burst (at 6 Gyr, decline times 0.1 Gyr). The spectra show the 
galaxies in their post-starburst phase at their maximum EW(H<S) and 1 Gyr 
after the burst. 



faster than in the Sd model. In contrast to the Sd model, the Sa 
galaxy at 6.4 Gyr still has an EW([On])> 5 A, as shown in Tab.[6] 
Its spectrum 1 Gyr after the burst looks like a typical passive spec- 
trum with a superposition of A-star features. The UV flux is very 
low and the Lyman-ff-line at 1216 A has disappeared. No emission 
lines can be seen, while the B aimer lines and the 4000 A break are 
present, even though EW(H(5) < 5 A for the Sa galaxy model (see 
Tab. [6). From the continuum ratios as well as from the strengths of 
the B aimer lines and the 4000 A break the different burst scenarios 
are distinguishable. 

The spectra in Fig.|13|show the galaxy models after the E+A 
phase, at a galaxy age of 9 Gyr. Now, 3 Gyr after the burst both 
spectra look like true passive galaxy spectra and the B aimer lines 
have disappeared. We conclude that studying E+A galaxies over 
a larger wavelength range, particularly including the UV range, 
promises a better discrimination of progenitor galaxies and trans- 
formation scenarios than only looking at EW(H6) and EW([On]). 
Our results also indicate that the analysis of Spectral Energy Distri- 
butions (SEDs) extending from UV/U trough optical and eventually 
NIR passbands will allow to encompass all post-starburst phases, 
even those, which are excluded from the classical E+A definition 
based on the optical region only. This will be investigated Paper II. 
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Figure 13. Spectra of an Sa model (burst 50%) and an Sd model (burst 70%) 
after their post-starburst phase, 3 Gyr after the bursts (at 6 Gyr, decline times 
0.1 Gyr). 



3.6 Comparison with Template Spectra 

In Figs.[l4]to[T6] different ages of the Sa galaxy model with a burst 
onset at 6 Gyr, a burst strength of 50% and a decline time of 0. 1 Gyr 
are sho wn together with the E+A galaxy template from lKennicuttl 
dl992al) . Fig.[l4]shows the Sa model 0.4 Gyr after the burst. At this 
age, the model spectrum is still too young to match the template. 
The flux between 4000 and 5000 A still is too high. When we look 
at Fig. [15] we can see that 0.7 Gyr after the burst the flux decreased 
enough to match the template very well. A little later, 1 Gyr after 
the burst, the flux between 4000 and 5000 A is already too low 
as compared to the template continuum. Fig. [T7] shows the same 
Sa galaxy model 7 Gyr after the burst, compared to Kennicut's E 
galaxy template. The model has passed the E+A phase and matches 
the template of an E galaxy very well. 

In Figs. [T5] and [T9] we compare the Sd model with a burst of 
70%, an onset at 6 Gyr and a decline time of 0. 1 Gyr with Kenni- 
cut's E+A and E templates, respectively. Shown are the Sd model 
spectra at the times when they best match. 

It is obvious that, in general, the model spectra match the tem- 
plates very well, however not quite as well as the Sa model spectra. 
This result is inline with investigations of NGC 3 92 1 in the liter- 
ature. It is proposed that NGC 3921 is a merger remnant of two 
disk galaxies, Sa and Sc, with a burst 0.7+/-0.3 Gyr ago, which is 
expected to evolve into a normal elliptical galaxy on a timescale of 

Our comparison of E+A galaxy spectra several Gyr after the 
starburst with template spectra of E and SO galaxies shows that 
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Figure 14. Sp ectra of an Sa mod el shown 0.4 Gyr after a burst and an E+A 
template from Kennicutt 1 1992a) (Sa, burst 50% at 6 Gyr, decline time 0.1 
Gyr). 
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Figure 15. Sp ectra of an Sa mod el shown 0.7 Gyr after a burst and an E+A 
template from Kennicutt 1 1992a) (burst 50% at 6 Gyr, decline time 0.1 Gyr). 
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Figure 16. S pectra of an Sa mo del shown 1 Gyr after a burst and an E+A 
template from Kennicutt 1 1992a;) (burst 50% at 6 Gyr, decline time 0.1 Gyr). 
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Figure 17. S pectra of an Sa m odel shown 7 Gyr after a burst and an E 
template from: K ennic utt 1 1992 a') (burst 50% at 6 Gyr, decline time 0.1 Gyr). 
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Figure 18. Sp ectra of an Sd mod el shown 1 .5 Gyr after a burst and an E+A 
template from Kennicutt 1 1 992a) (burst 70% at 6 Gyr, decline time 0. 1 Gyr). 
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Figure 19. S pectra of an Sd m odel shown 7 Gyr after a burst and an E 
template from' K ennic utt 1 1992 a)) (burst 70% at 6 Gyr, decline time 0.1 Gyr). 
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Figure 20. Time evolution of B-R for Sd models with different burst 
strengths (Sd, burst beginning at 6 Gyr, decline time 0.1 Gyr). 




Figure 21. EW(H(5) versus B-R of two Sd models with different decline 
times (Sd, burst 70% at 6 Gyr). 



E+A evol ve into early- t ype g alaxies in agreement with investiga- 
tions from iBicker et all <2003t) of the origin of SO galaxies in clus- 
ters. 



3.7 B-R Color of Post-Starburst Galaxies 

In Fig. [20] three Sd models with the same onset of the burst at 6 
Gyr, the same decline timescale of r=0. 1 Gyr but different burst 
strength (70%, 30% and 0%, i.e. SF truncation) are presented as 
well as an Sd model with an onset at 6 Gyr and a decline timescale 
of t=1.0 Gyr with a burst strength of 70%. For comparison an 
undisturbed Sd model is shown. We recall that the scenario with 
the long decline timescale does not lead to an E+A phase because 
it shows significant [On]-lines during its Hd-strong phase. A burst 
causes the galaxies to become very blue before they get red due to 
the truncation/termination of SF. The models with SF truncation or 
termination only quickly become red immediately after the onset. 
Therefore up to 12 Gyr the models with a starbust before the halt- 
ing of SF are bluer. Late type spiral galaxies are in general bluer 
than early type galaxies and therefore take longer to become red. 
For old ages all burst and SF truncation/termination scenarios con- 
verge to similar colors as the age differences between their stellar 
populations become more and more negligible and the photometry 
of the galaxies is dominated by stars of the Red Giant Branch and 
Asymptotic Giant Branch. 

The decline timescale of the burst has a strong effect on how 
fast the galaxy model becomes red after the burst. The model with a 
long decline time of 1 .0 Gyr is bluer than the undisturbed model for 
almost 2 Gyr, while the model with the shortest decline timescale 
of 0.1 Gyr is already redder than the undisturbed model after about 
0.5 Gyr. The symbols mark 1 Gyr intervals. In Fig.|2T]the EW(H<5) 
is plotted against the color B-R. With this plot it is possible to 
investigate the evolution in color during the Hd-strong phase. It is 
possible to see which galaxy models become red during their H<5- 
strong phase and which remain blue. Since blue as well as red E+A 
galaxies have been observed this inves tigation is of great interest 
dYang et al ] |2008l ; IZabludoff et alJI 19961) . The horizontal line is the 
5 A threshold, separating H<S-strong from non H<5-strong models, 
while the vertical line at B-R=l.l separates blue (to the left) from 
red galaxies (on the right). 

The value of 1.1 for the limit was chosen to lie between the 
values of the undisturbed Sb and Sc models at a galaxy age of 13 



Gyr, i.e. at redshift z=0. Model in the upper right square are red 
during their HrJ-strong phase. The arrow points into the direction in 
which the models evolve with time. 

In Fig.[2T]two Sd models with the same burst strength of 70% 
and an onset at 6 Gyr but different decline times r (0. 1 and 1 .0 Gyr) 
are shown together with the undisturbed Sd model. Again it has to 
be remarked that the model with the long decline time is not an 
E+A galaxy. The burst models separate from the undisturbed model 
after the 6 Gyr mark due to the rapid change in color at the onset of 
the burst. At the beginning of the burst, the color starts to change 
faster than the EW(Hr5). The model with the decline time of 1.0 
Gyr falls already under the 5 A threshold before it reaches the red 
side, i.e. it spends all of its H<5-strong phase in the blue. The model 
with a shorter decline time of 0. 1 Gyr crosses the blue-red border 
while still Hr5-strong. It is observable both in a blue and in a red 
H<S-strong phase because the value for EW(H<5) of the model with 
decline time 0.1 Gyr reaches a high peak value for EW(Hr5) and 
drops slowlier below 5 A than the color evolves towards the red. 
The two models in Fig. [2T] are representative for all models of our 
grid. Hence all E+A galaxy models with short decline times become 
red, while models with long decline times (not E+A!) always stay 
blue and do not cross the color threshold. 

For models with a decline time of 0.3 Gyr it depends on the 
time when the burst occurs wether they get red during their HrJ- 
strong phase or not. All models with a decline time of 0.3 Gyr and 
an onset at 1 1 Gyr become red during their HS-strong phase, while 
models with a decline time of 0.3 Gyr and an onset at 9 and 6 
Gyr can show either behavior. Some of these models have stronger 
values for EW(H<5) and become red during their H<5-strong phase, 
while others have EW(H<5) only marginally above the 5 A threshold 
and do not stay H<5-strong long enough to have a red H(5-strong 
phase. 

We conclude, that the later the type and the longer the decline 
time, the longer it takes for the models to become red. Red post- 
starburst galaxies are possibly a later stage of blue post-starburst 
galaxies with a strong burst, i.e. a merger or gas stripping sce- 
nario in a late stage. Scenarios that terminate the SFR on a long 
timescale, like those associated with harassment and strangulation, 
or scenarios with an intermediate decline time with a weak burst 
can be excluded for red H<S-strong galaxies. 

By looking at the colors of HS-strong galaxies we can con- 
strain their decline times. This result let us assume that these see- 
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narios are distinguishable by comparing the SEDs of post-starburst 
galaxies as investigated in Paper II. 



4 CONCLUSIONS AND OUTLOOK 

We investigated the transformation of galaxies when entering a 
high density environment like a galaxy cluster. In particular we 
focused on the E+A galaxy class, which contains galaxies in an 
important transition stage, a post-starburst phase, of galaxy trans- 
formation. E+A galaxies are defined to have strong B aimer absorp- 
tion lines but no significant emission lines. The spectrum of an E+A 
galaxy therefore looks like a superposition of the passive spectrum 
of an elliptical galaxy and deep B aimer absorption lines typical for 
A-stars. We examined the different processes of galaxy transfor- 
mation like mergers, harassment, gas stripping and strangulation in 
terms of their impact on the SF history of galaxies and identified 
the most likely progenitors and successors of E+A galaxies. For 
this purpose, we used the GALEV galaxy evolution code to calcu- 
late a grid of spiral galaxy models with different combinations of 
a starburst, i.e. a sudden increase of the star formation rate (SFR), 
and a halting of the SFR on different timescales, i.e. SF truncation 
and termination. With these scenarios it is possible to describe the 
effects on the SFR of the four different galaxy transformation sce- 
narios mentioned above. We varied the amount of gas consumed in 
the starburst (burst strength) and the galaxy age at which the star- 
burst, the SF truncation or termination begins (onset). The spectra 
and colors as well as the time evolution of Lick indices, in particular 
the EW(H<5) of the US absorption line at 4100 A and the EW([0 n]) 
of the [On] emission line at 3727 A as indicators of recent and 
ongoing SF were modelled and discussed. 

When investigating the time evolution of the EW(H(5), we 
found that the EW(H(S) is the result of a complex interplay of burst 
strength, decline timescale, onset of the burst and progenitor galaxy 
type: 

(i) The later the progenitor galaxy type, the easier an H(5-strong 
phase is reached. All Sd burst models become H(5-strong and the 
maximum value for EW(H<5) is reached by an Sd model. For a 
certain type of progenitor spiral galaxy, the essential parameter to 
cause an EW(H(5) above 5 A is the peak burst SFR i^bmst with dif- 
ferent values for each spiral type. 

(ii) The earlier the onset of the burst, the higher are the peak 
values reached by the EW(H(5). For bursts beginning at a galaxy 
age of 3 Gyr all scenarios are H(5-strong due to the larger available 
gas reservoir at such young ages. 

(iii) The longer the decline timescale of the SFR after a burst, 
the lower are the values reached for EW(H(5), but the longer is the 
duration of the H(5-strong phase. The US- strong phase can last up 
to 2 Gyr in cases of long decline timescales. 

(iv) The higher the burst strength, the higher are the peak values 
for the EW(Hcf). Only very few scenarios with a truncation or ter- 
mination of SF without a preceding burst have a Hc5-strong phase. 
Those are the late-type spirals with their relatively high SFRs and a 
truncation on a short timescale. They have EW(H(5) slightly above 

5 A for a short time. 

We find a significant effect of the metallicities in galaxies on 
their EW(H<5) and conclude that a higher number of E+A galax- 
ies at high redshift cannot be ascribed solely to a higher number 
of starburst galaxies but also to the lower metallicities of galaxies 
that cause longer Hc5-strong phases in the early universe and to the 
higher burst strengths enabled by the higher gas content. 



A striking result are the values for EW(H(5) that we found for 
undisturbed galaxies. We discovered that galaxies, which only fol- 
low the normal evolution, all go through a phase between ages of 
1.5 and 6 Gyr in which they have strong H(5-lines. Not taking this 
result into account in analyses of high redshift galaxies leads to 
an underestimation of emission lines, which are influenced by the 
Balmer absorption lines, for example the Ha-line and Hfi, and of 
the SFR derived from those. 

By investigating the second criterion for E+A galaxies, the 
[On] emission line at 3727 A, we found that the selection crite- 
ria for E+A galaxies commonly used exclude a significant number 
of post-starburst galaxies, in particular those galaxies with a long 
decline time. The conventional definition of an E+A galaxy is too 
narrow to encompass the full range of post-starburst galaxies. Pos- 
sible progenitors of E+A galaxies can only be galaxy mergers and 
maybe galaxies with gas stripping with a starburst, while the pro- 
cesses of harassment and strangulation can be excluded due to the 
long decline timescales they induce for SF. 

The comparison of E+A galaxy spectra several Gyr after the 
starburst with template spectra of E and SO galaxies showed that 
the successors of E+A galaxies are early-type galaxies. From our 
model spectra of E+A galaxies, we predict that studying E+A 
galaxies over a larger wavelength range, particularly including the 
UV, yields a better discrimination of progenitor galaxies and trans- 
formation scenarios. The analysis of Spectral Energy Distributions 
(SEDs) will allow to encompass all post-starburst phases, even 
those which are excluded from the classical E+A definition based 
on optical regions only. From our investigation of the time evo- 
lution of the B-R colors of our post-starburst galaxy models we 
found, that the progenitor galaxy type and the decline time have 
the greatest effects on the color evolution. The later the type and 
the longer the decline time, the longer it takes for the models to 
become red. Red E+A galaxies and red post-starburst galaxies in 
general are possibly a later stage of blue post-starburst galaxies 
with a strong burst and a short or intermediate decline time as e.g. 
triggered by a merger. Galaxy models with an intermediate decline 
time (>0.3 Gyr) with a weak burst and models with a long decline 
time (>1.0 Gyr) no matter how strong a burst was do not become 
red during their H(5-strong phase. Our results from examining the 
spectra and the colors of post-starburst galaxies indicate that it is 
possible to select and distinguish post-starburst galaxies by look- 
ing at their Spectral Energy Distributions (SEDs), which will be 
presented in Paper II. 

We have shown that E+A galaxies by definition include only 
transformation processes on a short timescale. The full range of 
post-starburst galaxies have not yet been investigated observation- 
ally. This needs to be done in order to achieve a comprehensive 
understanding of the galaxy transformation processes in high den- 
sity environments. 
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